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^r) . We present two-loop calculations on the axial- vector form factors Fa of semilep- 

tonic radiative kaon and pion decays in chiral perturbation theory. The relevant 
dimension-6 terms of the lagrangian are evaluated from the resonance contribution 
Ph. and the results of the irreducible two-loop graphs of the sunset topology are given in 

Q_i| detail. We also explicitly show that the divergent parts in Fa are cancelled exactly 

as required. 

> 

X ■ 1 Introduction 

Chiral perturbation theory (ChPT) [1, 2] has established itself as a powerful effective 
theory of low energy interactions. While it works for the strong interaction, it also in- 
cludes the electroweak one whose dynamics can be completely fixed by introducing the 
corresponding gauge bosons through the usual covariant derivative. Since the external 
momenta and quark masses are the expansion parameters for the generating function in 
ChPT [3, 4, 5], they need to be small compared to the physical scale of the chiral sym- 
metry breaking, i.e. about 1 GeV. Therefore, one expects that the semileptonic radiative 
kaon (pion) decays of K + — > / + z^7 (ti + — > / + z^7) can be well described in ChPT [6, 7]. It 
is known that these radiative decays [8, 9, 10, 11] could provide us with information on 
new physics [12, 13] by searching for the lepton polarization effects, which depend on the 
vector and axial- vector form factors F V a of the structure dependent parts. 

In this paper we deal mainly with the SU(3) ® SU(3) chiral symmetry. We will 
present two-loop calculations on the axial-vector form factors in K + — > {K^) and 

7r + — > Z + z/£7 (vri2 7 ) by virtue of the recent progresses in the p 6 -Lagrangian [14, 15, 16] and 
the massive two-loop integrals [17, 18]. Some remarks related to the form factors in 
and 7T; 27 follow: 

* Deceased 
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• The usual one-loop ChPT for the timelike form factor does not satisfy the unitarity 
or the final-state theorem and makes poor approximations [19, 20, 21]. Substantial 
corrections are expected at a higher order, especially at the two-loop level. 

• The form factors in the decays receive the first non- vanishing contributions at 0(p 4 ) 
but the first sizable ones, as well as the estimates of the accuracy, arise at 0(p e ). 

• The question of convergence in ChPT needs to be clearly addressed [22]. 

• In ChPT, the 0(p 6 ) contributions to the vector form factors Fy in the decays have 
been studied in Ref. [23], but those to the axial- vector ones Fa have been done only 
for Ti — > Ivf) based on the SU(2) x SU(2) symmetry [24]. 

In this paper, by using the relevant dimension-6 terms of the lagrangian [18, 20] from 
the resonance contribution [23, 24, 25], we perform a detailed calculation for the irreducible 
two- loop graphs of the sunset topology [21], which give the dominant contributions to Fa 
in both ti and K decays at 0(p 6 ). For completeness, we will also evaluate Fy and compare 
our results with those in Ref. [23] . 

The paper is organized as follows. In Sec. 2, we give the matrix elements for the 
decays. We review the lagrangians of ChPT to 0(p 6 ) in Sec. 3. In Sec. 4, we display 
the two-loop calculations for both vector and axial-vector form factors with the detailed 
formulas placed in Appendices A, B and C. In Sec. 5, we show the analytical results. We 
present our numerical values and conclusions in Sec. 6. 



2 The matrix elements 



We consider the decay of P(p) — > l + (l l )vi(s l ) r y(k) with P = K + or tt + , where 7 is a real 
photon with k 2 = 0. The matrix element M for the decay [6, 26, 27] can be written as 



M, 



' Gf 6 p M^( P , ky*(k)u(s')Y(l - 75MO , 



V2 



(1) 



where is the photon polarization and 9r+(-k+) — cos 9 (sin 9) with 9 being the Cabibbo 
angle. In Eq. (1), the hadronic part of the quantity M^ v is given by 



M^(p,k) = Jd 4 xe^(0\T(j; m (x)j: k (0))\p(p)) , 
which has the general structure 



(2) 



M^ u (p,k) = -V2Fj 



(P ~ k) u 



P(p-k) j; m P(p)) + V2F P g, u 



(p - kf - M 2 P 
-F A [(p - k)f,k u - g^k ■ (p - k)\ - r A (k ti k v - g^k 2 ) 
+iFve^ a pk a pP 



(3) 



where the first line represents the Born diagram, in which the photon couples to hadrons 
through the known KK'-f (7^7) coupling, with Fp being the P meson decay constant, and 
the subsequent lines correspond to axial-vector and vector portions of the weak currents 
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with Fy(A) being the vector (axial- vector) form factor. In Eq. (3), T\ is non-zero only 
for those processes with virtual photons, such as P~ — > l + l + l~Vi. In terms of the form 
factors Fa and Fy, we can write the vector and axial- vector parts in Eqs. (1) and (3) as 

M V (P -> l+ ua ) = e -9^ Fv£ ^* rk y^ (4) 

M A (P - l + m)=i—^-F A l v [k v {p-e)-g lu/ {p-k)ef 1 *], (5) 
respectively, where 

r = u(aV(l- 75 M0- (6) 

Both Fa and Fy are real functions for q 2 = (p — k) 2 below the physical threshold, which 
is the region of interest here, based on time-reversal invariance, and they are analytic 
functions of q 2 with cuts on the positive real axis. One of the reasons to perform the 
present calculation is that the q 2 dependence of the form factors starts at 0(p 6 ). 



3 The Lagrangians of chiral perturbation theory 

In the usual formulation of ChPT [1, 2, 28] with the chiral symmetry SU{3) L x SU(3) R , 
the pseudoscalar fields are collected in a unitary 3x3 matrix 



U{x) = exp 



.$(x) 



(7) 



where F absorbs the dimensional dependence of the fields and, in the chiral limit, is equal 
to the pion decay constant, F w = 92.4 MeV. The <E> is given by the 3x3 matrix 



( tt° + ^ V2n 





\ V2K- 



2F° 



V2K+ 



(8) 



where X, 



1, 2, • • • , 8) are the Gell-Mann matrices. 



An explicit breaking of the chiral symmetry is introduced via the mass matrix 

I ml 



X 





V o 












2m 2 K 



(9) 



m: 



where m n (K) is the unrenormalized 7r (K) mass. We note that the mass of rj to this order 
is given by the Gell-Mann-Okubo relation 



-m 



K 



1 2 



(10) 



The mass term in Eq. (9) is related to the quark masses by x — const ■ diag(m u , m^, m s ) 
with m u = 711,1. To calculate the form factors, we have to include the interaction with 
external boson fields. As previously stated, the electroweak gauge fields and are 
introduced via the covariant derivative 
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D„U = d^U + iUl^-ir^U, 

( 
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eA,Q ^_ 



cosOW- 
V sin^W- 



cos 6W+ sin 0W+ 








where Q is the quark charge matrix in units of e = g sin 9w, with 9w standing for the 
Weinberg angles and Gf/V% = fi ,2 /(8M^ / ). The final lepton pair (for clearity we will 
always refer to l + vi coming from W + ) appears in the leptonic current when substituting 



W, 



» I 



2V2M^ 



(12) 



The Lagrangians of ChPT contain both normal (or non-anomalous) and anomalous 
parts. Since the form factor Fy is related by an isospin rotation to the amplitude for 
7T° — > 77, it can be absolutely predicted from the axial anomaly. For this reason, we 
must also include the effect of the axial anomaly. At the two lowest orders, the full 
no n- anomalous Lagrangian is given by [1, 2] 



4 2) 

r(4) 



(13) 



+ L 2 Tr(D^U D u U ] )Tr(D»U D U U ] ) 



L x [Tr(D^UD»U r ) 

+L 4 Tr(D tl UD^)Tr( X U 1 + U X V ) 
+L,Tr{D ll UD»U\ X rf + U X *)) + L 6 [Tr( X U^ + U X r )f 

+L 7 [TrtfU - tf X )] 2 + L 8 Tr( X tf X tf + UxW) 
+iL 9 Tr(L llu D' 1 U D U U^ + R^rf D"U) + L w Tr(L^UR^) , 

where L^ v and R^ v are the field-strength tensors of external sources, defined by 



(14) 



J fJLU 



R 



(15) 



and {Li} are unrenormalized coupling constants. At 0(p 6 ), the non-anomalous Chiral 
Lagrangian contains 90 independent terms plus four contact terms for SU (3) [18] . The 
terms relevant to ijn-y) decays are found to be 



4 6) = yi7 (x + VO + yw (x+> ( MO + wi <x+ WD 

+2/82 (x+) (f+^f+) + ^2/83 (/+/«/ {x+ ) w^m^}) + iy SA (x+) (f+^u 1 

+W85 (f +IX uU t " X+ U U ) + it/ioo {f +l xv [f- P , /#] ) + 2/102 (X+f-^f-") 
+2/103 (X+) (/W^"> + ^104 (/+/- If^, X-]) + 2/109 (V^V'/H 



(16) 
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where 

u, = i{u\d^-ir^U -Uidp-itjtf} , 
X± = U^ + U^U, 

X±» = U^D, X U ] ±UD^U = V mX± -^{Xt>M • (17) 
The covariant derivative V M X = d^X + [r^, X] is defined in terms of the chiral condition 

IV = \ {U\d, - irjU - U(d„ - il^rf) . (18) 

The anomalous part begins at 0(p i ) with the Wess-Zumino (WZ) term C wz [29] 
containing pieces with zero, one and two gauge boson fields. The terms with one and two 
gauge bosons as well as the anomalous p 6 -Lagiangian [20, 29] are given by 

£$z,i = -^e^HriUd^dvUdaUkp-rfdpUdvrfdaUrp), (19) 

£$z,2 = --^e^ a Hr{d^dJ a Urp-d,Ud u r a UHp) 

+Ud^(£ u dJp + dJJ p ), (20) 

4«) = tC^ix-U.J+a^ + tCne^ix-lU^f-*?]) 

+C 22 s^ (^{V 7 / +7 „ f +aP }) + ■■■. (21) 

From the above expressions, the Feynman rules can be derived by expanding U = 
exp(iQ/F) everywhere in C = + + and identifying the relevant vertex mono- 
mials. In the next section we display the main result of the two-loop calculations for the 
form factors. 



4 The Form Factors 

To 0(p 6 ), the finite matrix elements in ChPT are obtained by multiplying the unrenor- 
malized Feynman diagrams obtained from C = + + with a factor y/Z per 
external meson, where Z is the wave function renormalization constant. To get these 
results, we start to calculate the mass and wave function renormalizations as well as that 
of the pion decay constant. 

Since the form factors at 0(p 4 ) are related to the finite counterterm contribution, we 
only need the wave function renormalization, m 2 Kn and F K n to 0{p A ) in our calculations. 
Explicitly, we have 



r-l 



1 - 



3F 2 
1 

4F2 



I(m 2 K ) + 21{ml) - 2A(2m 2 K + m')L 4 - 24m^L 5 



(22) 



I(m 2 ) + 2I(m 2 K ) + I {ml) - 32(2m 2 K + ml)L 4 - 32m 2 K L 5 , (23) 
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5m 7T = 



5rriK = 



Fir = 



where 



Q F 2 



-mll{m 2 v ) + 3mll(ml) - 48ml(2m 2 K + 
+96L 6 ml (2m 2 K + m 2 n ) + 96L 8 m 4 n 
7^ [4m 2 K I(m; 
+192L 6 m 2 K (2m 2 K + 



96L 4 m 2 K (2m 2 K + m%) - 96L 5 m 4 K + 192L 8 m 4 K 



F {1 + 
I U + 



1 
1 

+4m 2 K L 5 ]} , 



I{m 2 K ) - I (ml) + A(2m 2 K + ml)L 4 + 4m£Z, 5 



I(m 2 v ) 



\l{m 2 K ) - h{ml) + A{2m 2 K + )L 4 



A-D 



d D q 



m 



(2tt) d q 2 



m* 



16tt 2 



+ 1 + ln(47r) 
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In 



m 



(24) 

(25) 
(26) 

(27) 
(28) 



is the standard tadpole integral. It should be noted that the renormalization constants 
Sm(= rriphys — m) and 8F K;7r (= F K ^ — F) defined above are finite. The divergences and 
scale dependences of the loop integrals are cancelled by similar factors for Li in C^. 



4.1 The vector form factors Fy 

Using the chiral Lagrangians mentioned above, one immediately obtains the tree-level 
contribution for the anomalous parts of our processes. For the semileptonic radiative K 
decays, we have 



W,K,tree 



4V2n 2 F 



1 _ 256^^ + 2^2^ _ m 2 )C W + ^2^2 + k 2 )C 



64 



(29) 

Loop corrections to the above tree-level contribution proceed through diagrams involving 
at least one vertex given by the WZ lagrangian. As shown in Figure 1 with P = K, due 



K-( P ) 




Figure 1: One- loop diagrams that contribute to Fy in P^ 7 



to the initial order, most one- loop diagrams can contribute to Fy to 0(p 6 ) constructed 
with one vertex coming from the WZ Lagrangian and another one from the lowest order 
chiral lagrangian. Performing the calculation in the three-flavor case, we get 

FvKloov(a) = 7= ' 7= M m l + 2?77/L + 3m 2 ) 

V,K,loop(a) 4 ^ n 2 F l 67r 2 (v /2 F )2 L V K «> 
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m; 



- 3 < ln h# 



-4m 2 + ^ 2 )A 



V,K,loo P ( b ) 4v ^ 7r 2 F • 167r 2 ( ^2 F) 2 

+ 4: J ( m l ~ ^(1 — 



In 



m 2 - x(l - 



U2 ~ 



V,K,ioop(c) 



4-\/27r 2 F 16t: 2 {V2F) 2 



'-ml -m\- 2m 2 K 



' xm 2 + (1 — x)m\ — x(l — x)q 2 



+2 J (xm 2 + (1 — x)m 2 K — x(l — x)q 2 ) In 

r. /*/ 2 /-, x 2 /-, \ 2\ , Z^™* 2 + (1 - x)m 2 K - x(l - x)q 2 \ 
+2 J (xml + (1 - x)m\ - x(l - x)q 2 ) In I £ ^ 1 CL. j dx 



dx 



from the three loops in Figure 1, respectively, where 

A = i + 1 + ln(47r) - 7 . 



, (30) 



(31) 



Putting all of the contributions together as well as the usual renormalizations of the 
pseudoscalar wave-functions and decay constants, we obtain the following results: 



V,K 



L— {l - ^n 2 m 2 K C 7 + 2567r 2 (mK - ml)C n + ^vr 2 (g 2 + k 2 )C 22 
n 2 F K y 3 3 



4^/271 ~r k 
1 



5 - ml)\ + \{k 2 + q 2 )\ - ^m 2 In ( ^ ) - 3m 2 , In ( ^ 



ml 



16n 2 {V2F K ) 2 
+2 y xm 2 + (1 — x)m 2 K — x(l — x)q 2 



3 V " ^ 3 

m 2 — x(l — x)k 2 



In 



+2 



xm 2 + (1 — x)m 2 K — x(l — x)q z 



In 



In 



ml — x(l — x)k 2 \ 
xml + (1 — x)m 2 K — x(l — x)q 2 
'xm 2 + (1 — x)m|- — x(l — x)g 2 



/// 



Similarly, for ix^ we derive 
1 



Fv ' n AV2n 2 F n 



256 o 2^y 64 2/ , ^.^ 1 



2 2 

-4m 2 In - 4m 2 In ^ + 4 / (m 2 , - x(l - x)k 2 ) In 

+4 / (m 2 - x(l - x)q 2 ) In < ~ < l ~ x )<? + h k 2 + q 2 )\ 
J 3 



(32) 



(33) 



We note that the numerical results of Eqs. (32) and (33) agree with those in Ref. [23] as 
will be shown in Sec. 6. The one-loop corrections contain divergent terms proportional 
to A = - + 1 + ln(47r) — 7 coming from the dimensional regularization scheme. Obviously, 
the presence of these divergent terms requires the introduction of the corresponding coun- 
terterms in the anomalous section of the Lagrangian at 0(p 6 ), which were already given 
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in Refs. [20, 30, 31]. Their infinite parts cancel the A-terms in Eqs. (32) and (33) and the 
coefficients are simply substituted for by the remaining finite parts, the renormalized 
coefficients Cf r . The values of these finite contributions from the counterterms to our 
processes are not pinned down in ChPT and they have to be deduced from data fitting [29] 
or, alternatively, from the hypothesis of resonance saturation (RS) of the counterterms 
[23]. The estimations of using RS are described in Sec. 6. 



4.2 The axial-vector form factors Fa 

In this subsection we aim at the extraction of the form factor Fa, which is the only one that 
has a contribution proportional to g^ v {p • k). The presence of g^ v requires that the axial- 
vector and vector insertions are in the same one-particle irreducible subdiagrams. This 
immediately removes a large part of the diagrams. Furthermore, the (p ■ k) kinematical 
factor guarantees that it is not part of the internal Bremsstrahlung contribution. We now 
discuss the contributions from the diagrams to 0(p 6 ). 



4.2.1 Tree-level diagrams 

With the chiral Lagrangians in Eqs. (13), (14) and (16), one obtains tree-level contribu- 
tions for the processes as follows: 



4j\/2, 
~F 



F A%,tree ~ —^r( L 9 + ^lo) , (34) 



F A % ree = ^(L 9 + L w ), (35) 



(6) 48ygm| 48y/2(2m* K + mp \^m\ 

^A,n,tree ~ ~Vn p 2/18 p r 2/81 ^ 



\6V2{2m 2 K + ml) \^m\ 8^/2(2m\ + m 



,2' 



+2/82 J, 2/83 J, 2/84 

— 1/85 p r 2/100 p 2/102 p 

16V2(2m 2 K + ml) , 16y/2ml , 8^(2^ • k - pi ) 

-2/103 p 1" 2/104 p \- 2/109 p 

-2/110 p , (36) 

A8V^m 2 K 48V2(2m 2 K + ml) \§^2{2m 2 K + ml) 



F A,K,tree = ~2/l7 ^ ~ ~ 2/18 ~ + 2/81 

16V2(2m 2 K + mi) 16V2(2m 2 K 

+ m i) 

+2/82 p 2/83 7^p 

8V2(2m 2 K + ml) 8^2{Am 2 K - ml) 

_ 2/84 p 2/85 ^p 

8y/2{p 2 K -p K -k) 16V2m 2 K \^(2m 2 K + m 2 n ) 
+2/100 p 2/102 p 2/103 p 

16V2m 2 K 8V2(2p K -k-p 2 K ) A^2p K ■ k 
+2/104 =; + 2/109 f 2/no f • ( 37 ) 
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We note that for the tree- level contributions in Eqs. (36) and (37) at 0(p 6 ) one needs to 
perform renormalization with the finite parts, which will be discussed in Sec. 5. 



4.2.2 One-loop diagrams 

We now consider the one- loop diagrams shown in Figure 1. Since it contains only photon- 
even- meson vertices in the non-anomalous chiral lagrangian to 0(p 4 ), Figure 1(c) does 
not contribute to Fa- Moreover, one-loop diagrams with an 0(p A ) vertex insertion on a 
propagator in the loop never produce the factor (p- k) and hence do not contribute to Fa- 
From Figures 1(a) and 1(b), we get 



FA,ir,loop(a) 



1 A,K ,loop(a) 



F A ,ir,loop(b) 



FA,K,loop(b) 



3F 3 
2Lio 



\uV2I(m 2 K ) + 28v / 2/( m„)\ 
l ,., [l0V2I(m 2 K ) + 20V2I( m„)\ , 

[6V2I(m 2 v ) + 2AV2I(m 2 K ) + 12 v / 2/( 



3F 3 
2Lio 



3F 3 
-16v^^i 
F 3 



[3v^/(mJ) + \%y/2I(m 2 K ) + 9V2I(ml) 



Kml) + 



I(m 2 K ) + 2I(ml) 



Kmi) - 



F 3 

IQV2U 2 U2U 2 , 



F 3 
2\/2L 9 



F 3 



F 3 

I(ml) + 2I(m 2 K ) 



8V2L 3 
F 3 



8V2L 3 
F 3 



I (ml 



1 

2" 



J(m^) + -I (ml) 



(38) 



4.2.3 Two-loop diagrams 

The two-loop diagrams which may contribute to Fa are shown in Figure 2. The last six 
diagrams with nonoverlapping loops in Figure 2, which can be written as the products 
of one-loop integrals and produce no (p ■ k) factor, do not contribute to Fa- The only 
possible non-vanishing diagrams are the first three irreducible ones in Figure 2. 

Since there are three different mass scales of (m^,m K ,m r ^) with the same order of 
magnitude in the SU(3) ChPT, the irreducible integrals can no longer be expressed by 
elementary analytical functions. We will quote only the numerical results in Sec. 5 and 
Sec. 6. We now give the detailed calculations for the g^ v terms of Ma in Eq. (5). It 
is clear that the first irreducible diagram in Figure 2 does not contribute to Fa since 
there is no (p ■ k) term. The second and third irreducible diagrams with genuine massive 
two- loop integrals [21] are depicted in Figures 3 and 4, respectively. Our calculations on 
these diagrams are summarized in Appendix A. 

We point out that the two-point irreducible diagrams in Figure 3 vanished in Ref. [24] 
for 7r^ 27 in the linear sigma model parametrization in which there is neither a photon nor 
a four-pion vertex [32]. As shown in Sec. 5, these diagrams play a very important role in 
the cancellation of divergent parts with 1/e 2 in the final results. 

We note that the two-loop amplitudes shown in Appendix A have been classified by 
the functions {/, A, B...L} related to each diagram. For clarity, we always refer to all 
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Figure 2: Two-loop diagrams for the form factors in Pn-y- 



mudf, d„, 1711,8^ mi, 9^ 




m 3 m 3 mz 

I(mi,m 2 ,m 3 ) II(mi,m 2 ,m 3 ) III(mi,m 2 ,m 3 ) 



Figure 3: Two-point irreducible diagrams 

inward particles for vertices. The Feynman diagrams which contain tensor structures at 
the numerator can be reduced to the calculation of scalar integrals P^\ a2 a3 by introducing 
the transverse components of the loop momenta s and q as s± = s M — (s ■ l/l 2 )l^ and 
q± = q 11 — (q ■ l/l 2 )^, respectively, for the two-point diagrams [21]. For the three-point 
diagrams, we first combine the denominators by using the Feynman formula 1/ab = 
Jq dx/[ax + 6(1 - x)f and then replace I by V = I + kx, i.e., = s>* - (s ■ I'/l' 2 )^ and 
Q± = q^ — iq-l'/l' 2 )^- The corresponding relations among {/, A, B...L} and P^\ a2 a3 
are displayed in Appendices B and C , where the scalar integrals i^ a 2 a 3, defined in the 
Euclidian space, are given by 

P:U„Mu m,, m.; P) = / ^ ( ^ + mj) . 1(g3 ( ;^ + 8)a + mi) ., ■ (39) 

The detailed contributions from all diagrams for Fa in n^-y and are summarized in 
the next section. 
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m 3 

^(7711,7712,7713) 



m 3 

5(7711,7712,7713) 






J (mi, m.2, 1713) K(mi,m2,m^) L(mi,m2,mz) 

Figure 4: Three-point irreducible diagrams 

5 Analytical Results 
5.1 Renormalization scheme 

In our calculations, we use the following dimensional regularization and renormalization 
scheme [1, 2, 22, 28]. Each diagram of 0(p 2n ) is multiplied by a factor (cyu)( n_1 ^ D_4 ), 
where D = 4 — 2e is the dimension of space-time and c is given by 

Inc = -\ [1 - 7 + 1h(4tt)] - | + ^ + 0(e 2 ) . (40) 

Using the renormalization factor of (c/x)^ -1 ^- 0-4 ), the low energy constants Lj of C$ in 
Eq. (14) are defined by 
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Li = (c//) (D_4) £-;(//,£>) . 
Similarly, for the Cffl parameters yi in Eq. (16), 

yi = (c^ D -V yi (t,,D). 



(41) 



(42) 



We note that L^fx, D) have the same /i-dependences as the one-loop integrals, whereas 
yi(fi, D) behave like the two- loop ones. Their values at the two different scales of //i and fi 2 
are related by D) = (/i 2 / /i^ '^ Li(ii 2 , D) and y^fx^D) = (/i 2 / 'A t i) 2( ' D_4) ?/i(A t 2, D), 

respectively. 



5.2 Analytical forms of Fa 

We now try to obtain the analytical forms of F A from each diagram to 0(p 6 ) at the scale 
of m p . From Sec. 5.1, for the unrenormalized coefficients in the chiral lagrangian we 
use [22, 28] 



(D-4) 



li 



= V 



(D-4) 



32tt 2 ( 

li 



327T 2 



7s 



In 



m 



+ L r i (m p ) + eL\ £ \m p ,n) 



(43) 



which is a Laurent series expanded around e = 0. Of the values in Eq. (43), 7$ are shown 
in Table 1 [22], 7^ = — 1 — ln(47r) + 7, L\ correspond to measurable low energy constants 



and L\ £) are given by L\ £> (m p ,n) = -L[(m p )(7 s + ln(m 2 / > 2 )) + 7* ■ f(m p ,fj). Similarly, 



Table 1: Coefficients of 7$ in the Minkowski space [22]. 



i 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


li 


3 

32 


3 

16 





1 

8 


3 

8 


11 

144 





5 

48 


1 
4 


-1 
4 



for the chiral lagrangian Cffl, we have 



Vi = 



2(D-4) 



F 2 



+ 



rf + iTV) 



1 



,(2) 



327r 2 e 



/i 



2(D-4) 



F 2 



& r (™p) + 



(if> + r 



S L) («v 



32vr 2 



1 I m n 

--2 7s -21n 
e \ /i 



1 32 2 vr 4 e 2 

2 
p 
2 



(2) ^(m p ,/i,ej . 
1 32 2 vr 4 e 2 ' ' K ' 



the relevant constants of vf\ vf^ and vf^ are shown in Table 2 [22]. We note that 
f(m p , ji) and g(m p , /i, e), which receive contributions from the (e/2) term in Eq. (40), will 
not contribute to F A - 

By writing the unrenormalized contributions to Fa of 0(p 4 ) and 0(p 6 ) diagrams as 
F A ,tree(p 4 ), F A ,tree(p & ), F A1 _ loop (p 6 ) and F A , 2 -ioop(p 6 ) , respectively, using the wave function 
renormalizations, we find 



F A = 



F A ,tree(p 4 ) (1 + SF P ) + F A)tree {p & ) + F A ,l-loop{P & ) + F A2 _ loop (p e ) , (45) 
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Table 2: Coefficients of T- ' and T\ ' ' with the double-pole and single-pole divergences, 
respectively, in the Minkowski space [22] . 



Vi 


i 


i67r 2 rf } 


r(i) 

i 


17 


19 

64 


13 

768 


2 rr _j_ 4 rr _j_ 8 rr _j_ 3x r 
3-^1 + 3^2 + 9^3 + 4^5 


18 


67 

192 


1 

2304 


8rr_i 8 rr i 23 r r i 3 x r , lrr 
3-^1 T 9-^2 + 27-^3 + 2^4 + 4^5 


81 








3 r r A j r 
4-^9 4^10 


82 








lrr 1 T r 
4-^9 4^10 


83 


1 

4 


5 

384 


lrr -lrr 5 rr 3rr i 3 rr 
3-^1 6^2 3f> 3 4 5 "1~ 8 9 


84 


4 

— f- 


5 
144 


32 r r 32 r r 157 rr f?rr rr x I rr 
— — iv 1 — -g-^2 — lA^Z ~ 0iv 4 5 ' 2^9 


85 


2 


5 
192 


2 rr 11 rr 16 rr 3 rr i 3 rr 
3^1 — T^2 ~ "Tf-^3 — 2^5 4^9 


100 


5 

24 


49 
1152 


3 1 6 2 ~ 4 3 3 4 4 5 8 9 


102 


25 
64 


17 


2 r r 4 r r 8 r r 3 t r , 3 r r 
_ 3^1 ~~ 3^2 ~~ 9^3 ~~ 4^5 + 4^10 


103 


73 

192 


— m— 


8 rr 8 r r 23 rr 3 rr lrr i lrr 
3-^1 9 2 27 3 2 4 4 5 4 -^10 


104 


7 

96 


|3D4 
144 


irrilrr l rr _Lll r _Ll rr _LO rr 
~~ f> 1 + 12^2 ~~ 8-^3 "+" 6^4 "+" 8^5 "+" 16^9 


109 




3 
128 


lrr 
"2^9 


110 


6 


31 

576 


2 rr i lrr lrr i 2 x r i lrr i lrr 
3-^1 3^2 ~~ 2^3 "+" 3^4 "+" 2^5 "+" 4^9 



which leads to 
F A 

where P = n or K. In Eq. (46), for P = it, we have 



F A ,tree(p 4 )(l + SF P + ^5Z P ) + F Attree (p 6 ) + F AA _ loop (p 6 ) + F A ^ loop (p 6 ) ,(46) 



1 4\/2 
F Atree (l + 5F n + -5Z n ) = (Lg + L 10 ) 



1 - 



' I{m 2 



F 2 



4^2 



+ ^o) 1 



16tt 2 F2 



i-2 7 »-21„(^ 



A' 



1 



16vr 2 F 2 



3 \ m 2 



2m| ^ / mj 



mt 



2m 2 



(47) 



loop 



GV2F* 



% 2 



2% 




[6(2LJ - I/ 2 )m 2 + (m 2 K + 2m 2 )(3^ + 5Lg + bU w ) 



6V2F^n 



m K ln 




-12L^ + 6L^)m 2 ln 



ml 



- {3L r 3 + 5L r 9 + 5L r w )x 



+ 2m 2 ln 



m; 



(48) 



F A . 



2— loop 



6F3(2tt) 8 



7r 4 (116m 2 K + 184m 2 + 21m 2 ) h 
12y/2 e 



2 7s 
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3tt 4 



2V2 



^-2 7s -21n(-| 



(p-k) 



-1421.4 - 1167.7(p • fc) + 1228.0 + 1123.2(p • fc)} , 



(49) 



,iree 



e 



/ m„ 

21n U 



3tt 4 



6F3(2tt) 8 2v^ 



7r 4 (432m| + 531m 2 ) 



+ 



6F3(2tt) 8 36v/2 Gv^i^vr 2 

+6m 2 (2^ - L\ + LI + 2L r 9 + 2I7 W )] } 
4\/2 



3m^(L5 + 2LJ + 2LJ ) 



F 3 



{4m 2 , (6y r 18 - 2y r 82 + y r 84 + 2y r W3 ) 



+2m 2 (6y r 17 + 6y[ 8 - 2y r 81 - 2y£ 2 + 2y r 83 + y r 84 + y r 85 - y r 100 + 2y[ 02 
+%Io3 - 2 Vw4 + Vim) + pM 2 3/ioo - 4y[ 0g + y[ 10 )} . 
For P = K, we get 

1 /^K) 



1 4\/2 
F A , tree (l + 5F K + -5Z K ) = — ^— (Lg + L; o ) 



1 - 



F 2 



4 



/(mp I(m 2 ) 



--2 7s 



2ln U 



2\ 1 



X ,4 , < 

.4 2 4. 



16tT 2 F 2 



4 \ m 2 



+ 



2 1 m 2 



H -In 



(50) 



(51) 



loop 



4V2F*ir 2 



1 

--2 7s 
e 



/ m„ 

21n u 



{m^(8L^ - 4L r 2 + 4L r 3 + + 6L r w ) 



1 



4-\/2F; 



3 ^2 



A' 



7T" 



-(2L^ + 3L5 + 3^ )m 2 ln 



-(8L\ - AL r 2 + 4L^ + 6U 9 + 6L r 10 )m 2 K In 



m 



K 



- {Ll + L\ )ml In 



m, 
m 



(52) 



Fa,2- 



loop 



7r 4 (99mi- + 20m 2 - 12m 2 ) 



3tt 4 



4^ 



--2 7s 
e 



2^ 



1 / m 

--2 7s -21n [—: 
e \ fl- 



ip -k) 



-1465.2 - 923.0(p • k) + 1266.9 + 1208.3(p • k)} , 



/ m 2x n 
2111 / 



(53) 



.tree 



1 



- 2 7s - 2 In 



m 



1 



7r 4 (747m 2 f 



3tt 4 



6F|(27r) 8 2x/2 



(p- fc) 



216m 2 ) 



6^(2^ 



36^ 



+ 



AV2F^tt 2 



ml(2L r 3 + 3I7 9 + 3L[ ) 
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+m 2 K (8L r 1 - 4L r 2 + AL r 3 + 9U 9 + 9L r w )] } 

+2m^(18y[ 7 + 36y[ 8 - Ayl, - 12*4 + 4^ 3 + 6y r M + 4y r 85 - 3y r wo 
+6y r 102 + 12y[ 03 - Qy[ 04 + 3y r 109 ) + 3pk(2y[ 00 - 4y[ 09 + y[ 10 )} . (54) 

We note that, in the above expressions, one special case must be treated separately. 
For the finite part of F^-ioop, the functions of g(y) and fi(y) in Appendix B seem to 
introduce additional singularities in the integrals, such as at y = and 1. However, this 
problem can be resolved by noticing that at y = (the situation is the same at y=l), 
the function, e.g. g(y), behaves like m(a; 2 ), which is integrable. Consequently, the main 
question of evaluating the finite part of Fa,2-Ioo P is how to implement the formula in a 
computer program by correct and numerically stable forms. As seen from the last two 
terms for F^^-ioop in Eqs. (49) and (53), we have found the reliable and stable numerical 
results due to the contributions of functions hi, i.e., the integrals of g and f iy by fitting 
the numerical data with recursive analytic methods. 

We remark that in Eqs. (47)- (54) we have explicitly shown the single poles, subtracted 
via FA,tree{p 6 )- We emphasize that in our results there are no divergent parts with 1/e 2 
and all the terms related to 1/e are cancelled explicitly by the renormalization of the 
coupling constants in C$ as well as the Gell-Mann-Okubo relation in Eq. (10). It is 
clear that the disappearance of 1/e 2 terms relies on the two-point irreducible diagrams in 
Figure 3. Moreover, our results are scale independent since the scale terms with In ji 2 can 
be grouped into the ones with 1/e, i.e., they are always associated with 1/e terms. These 
also serve as checks of our calculations. Now we can extract the axial-vector form factor 
by placing the related physical quantities into the Eq. (46). Explicitly, we may write Eq. 
(46) into more transparent forms 



V = {66.86(LS + LSo)} 

+{(2.41 - 122.96^))% - 4.82% - 125.35% + 4.82% 

pi pi pi pi 

+ (-2.41 + 245.95(pA;))^f - 61.49(pfc)^ - 14.46^ - 376.05^ 

+4 . 82 |i + 125.35^ - 4.82|§ - 62.67^ - 2.4l|§} 

+ {12.30^ - 6.15L?, + 16.11LJ + 26.85L£ + 26.85L£ } 

+ {-1.70- 1(T 2 - 3.92- 10~ 3 (pA;)} (55) 



and 



F AK = {54.99(25 + LJ )} 



+{(24.75 - 101.01(p£;))% - 49.50% - 102.97% + 49.50% 

pi pi pi pi 

f V V T 

+ (-24.75 + 202.03(pA;))^f - 50.50(p£;)^f - 148.50^| - 308.9o|g 
+34.32^| + 102.97^| - 34.32^| - 51.48^| - 32.34^|} 
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+ {22.08L^ - 11.04L5 + 12.75L5 + 21.71LJ + 21.71LJ } 
+ {-0.97 • 10~ 2 + 13.93 • 10- 3 (j»A;)} , 



(56) 



where the four {• ■ •} terms in Eqs. (55) and (56) correspond to those in Eq. (46), 
respectively. 

6 Numerical values and conclusions 

As shown in section 5, the divergent terms for Fa in loop-diagrams are cancelled by the 
corresponding counterterms in the Lagrangian at 0(p 6 ). The infinite parts cancel each 
other and thus they can be simply substituted by the remaining finite part of the coun- 
terterms, y\. We now study the finite parts which contain the actual physical information. 
We will present the results in numerical forms, with the scale at m p = 0.77 GeV. In Table 
3, we show the standard values for the couplings L\ in C$ [6] and those in the two- loop 
calculation of ChPT; we chose the Main Fit in Ref. [33] as an illustration. Other two-loop 
studies in ChPT can be found in Refs. [34, 35]. We note that in the table the central 
value of a[ = —5.5 is kept, and our numerical results for Fa are sensitive to this value. 
Several 0(p 6 ) low-energy constants of the normal chiral Lagrangian have been evaluated 
from the RS. In Table 4 we illustrate the values of yi in the lowest meson dominance 
(LMD) approximation [25] and the resonance Lagrangian (RL) [23, 25, 36]. 



Table 3: Values of L\ in 



10 3 L[ 


1 


2 


3 


9 


10 


0(p 4 ) [6] 


0.4 ±0.3 


1.35 ±0.30 


-3.5 ± 1.1 


6.9 ±0.7 


-5.5 ±0.7 


Main Fit [33] 


0.53 ±0.25 


0.71 ±0.27 


-2.72 ± 1.12 


6.9 ±0.7 


-5.5 ±0.7 



Table 4: Values of y t in £(?). 



Ui(in units of 10 4 /F 2 ) 


2/100 


2/104 


2/109 


2/no 


LMD 


1.09 


-0.36 


0.40 


-0.52 


RL I 


1.09 


-0.29 


0.47 


-0.16 


RL II 


1.49 


-0.39 


0.65 


-0.14 



To study the vector form factors, we need to consider the anomalous chiral Lagrangian. 
The set of anomalous coefficients is treated by phenomenological fitting in ChPT as well 
as by the two main alternative models of vector meson dominance (VMD) method and 
constituent chiral quark model (CQM) [29]. The relevant terms for our purposes are shown 
in Table 5. Other physical inputs are itik = 0.495, m n = 0.14, m v = 0.55, Fk = 0.112, 
F = 0.0871 and F n = 0.092 GeV. We note that some of the actual values of F n and the 
masses in our calculations may be different from those in the literature. In particular F n 
could differ around 1% from one paper to another; however, we expect that the changes 
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Table 5: Values of Cf r in £(?) in various models [29]. 



C^ r [10- 6 GeV- 2 } 


7 


11 


22 


ChPT 


0.013 ±1.17 
20.3 ± 18.7 


-6.37 ±4.54 


6.52 ±0.78 
5.07 ±0.71 


VMD 






M ja 2 ^ 8.01 


CQM 


0.51 ±0.06 


-0.00143 ±0.03 


3.94 ±0.43 



on our 0(p & ) results due to the different sets of parameters are less than 5% since 0(p 6 ) 
contributions are at least proportional to F%. 

To compare our results with those in the literature, we use dimensionless form factors 
of fv,Ai defined by 

fi = mpFi, (i = V,A) (57) 

to replace Fy t A- 

In Figures 5 and 6, we plot the dimensionless vector and axial- vector form factors fv,A 
as functions of q 2 with the photon on mass-shell for 7r e27 and K e2l , respectively. Similar 
figures can also be drawn for the \x modes. In Table 6, we show the form factors of Ja 
at q 2 = at 0(p i ) and 0(p 6 ) with SU(2) and SU(3) symmetries as well as experimental 
values for P = K and it. 



Table 6: f A at q 2 = for P = K and tt. 



f A (q 2 = 0) 


0(P 4 ) [6] 


0( P e )\su(2) 


0(P 6 )\SU(3) 


Experiment 


P = K 


0.041 




0.034 


0.035 ± 0.020 [37, 38] 


P = n 


0.0102 


0.0117 [24] 


0.0112 


0.0116 ±0.0016 [39] 




0.005 0.01 
q 2 G eV 2 



0.015 0.02 



0.11 
0.1 
0.09 
0.08 




0.05 0.1 0.15 0.2 0.25 
q 2 G eV 2 



Figure 5: The vector form factors fy as functions of the momentum transfer q 2 for ix^ 
and K(2~t with I = e. The dot, solid and dashed curves stand for the contributions of 
OG» 4 ), 0(p 4 ) + 0(p 6 ) in VMD and 0(p 4 ) + 0(p 6 ) in CQM, respectively. 
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0.016 
0.015 
0.014 
0.013 
< 0.012 
0.011 
0.01 
0.009 



0.005 0.01 
q 2 G eV 2 



0.015 0.02 



0.045 
0.0425 
0.04 
0.0375 
0.035 
0.0325 
0.03 
0.0275 



0.05 



Kl2y 



0.1 0.15 

q 2 Ge V 2 



0.2 0.25 



Figure 6: Same as Figure 5 but for the axial-vector form factors j \. The dashed and 
solid curves represent the contributions at 0(p 4 ) in Table 3 and the fitting 0(p 4 ) + 0(p e ), 
respectively. 

In Figure 5, the dot, solid and dashed curves stand for the contributions to fy at 
O(jo 4 ), 0{p A ) + 0(p 6 ) in VMD and 0{p A ) + 0(p 6 ) in CQM, respectively. We note that 
for K[ 2 ~, in Figure 5 F has been set to F n in Eq. (29) for the curve of 0{p A ) as in the 
literature and F K in Eq. (32) for those of 0(p i ) + 0(p & ). As shown in Figure 5, the 
0(p 6 ) contribution obtained for the 7r radiative decays is very small (< 5%) for all the 
kinematical allowed values. However, it is interesting to see that the 0(p 6 ) correction for 
K^j is much larger. 

In Figure 6, for /a, the dashed and solid curves represent the contributions at 0(p 4 ) 
in Table 3 plus the two-loop calculations using either the LMD or the RL determinations, 
respectively. It is easy to see that, as shown in the figures, the two-loop contributions to 
] a are sizable and destructive compared with those at the pure 0(p 4 ) for both 7r and K 
modes, but their g 2 -dependences, which are dominated by the irreducible diagrams, are 
small. At q 2 = 0, the contributions to ] a from f'A,tree{p 6 ) are vanishingly small, which 
implies that the final results of Ja are insensitive to the known values of y[. However, 
those from the irreducible two-loop diagrams and the one-loop diagrams with one vertex 
of give the dominant corrections to Ja at 0(p 4 ). All together the 0(p 6 ) corrections 
keep around 25% for both decays. We remark that the uncertainties in Eqs. (55) and (56) 
due to the errors of y\ [25] are less than 1%. Moreover, our results of tree contributions 
at 0(p 6 ) are also insensitive to the choice of the scale, as expected. 

We note that, as shown in Figure 5, the numerical result for n — > eu e j at 0(p 6 ), using 
the SU(3)<S>SU(3) chiral symmetry, is found to be comparable to the one in Ref. [24] from 
resonance estimates of 0(p 6 ) low-energy constants based on SU(2)(B)SU(2). Furthermore, 
our result of the 0{p & ) correction for Ja in Ki 2l also confirms the speculation in Ref. [24]. 

In summary, we have studied the 0(p 6 ) corrections to the vector and axial- vector form 
factors in ir^-y and K^-y decays. These include the contributions from loop diagrams and 
the ones from higher dimension terms in the lagrangians. The former can be exactly 
calculated in terms of the known parameters of the chiral lagrangians. The latter is 
mainly evaluated from the resonance contributions. For the axial-vector form factors of 
fA, we have found that the divergent parts cancel order by order in ChPT, while the 
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finite numerical results in both K and it modes contain considerable corrections from 
loops diagrams; they also agree with the recent experimental determination [37, 38]. This 
demonstrates numerically the statement about the final-state theorem mentioned in Ref. 
[19, 20, 21]. Finally, we remark that our result of f a at q 2 =0 for the kaon case is consistent 
with that found in the light front QCD model [40, 41] as well. 
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Appendix A 

From Figure 3, we have 

M A ,n,2- P oint = — - 24i?3 cos 0[ -40V2I(m K , m K , m*) + 40V2I(m 7r ,m K ,m K ) 

+24V2I(m K , rriK, in v ) — 24\f2I(m r] , uik, tuk) — 240v / 2//(mx, rn^, uik) 

—40V2II(m w , uik, tuk) — 320\[21 1(171^, m n , m^) 

-48V2II(m K , m K , m v ) — 72V2II(m v , m K , m K ) 

+240V2I II(m K , m K , m^) + 40V2III(m K , m n , m K ) 

+320v / 2//J(m 7r , m w , m n ) + 96v / 2///(m^, mx, tuk) 

+72V2III(m K , m v , m K ) - 48V2III(m K , m K , m v )\, (58) 

M AKt 2-point = ^ 6 ^24F 3 sin 6l [ _36v/ 2 / ( m ^' "V m K) + 3§\/2I(m K , m v , m v ) 

3§\[2I(m K , 
- 336^21 1(m K 
-3QV2II(m K ,m n ,m n ) — lQ8V2II(m 7T ,in n ,mK) 
—36\/2II(m v , itik, m^) - 9QV2II(m K , m n , m„) 
— 12 \/2 II (m^, 171k, m v ) + 3Q\f2III(m ri , m^, m^) 
+336\/2III(m K , m K , m K ) + lQ8V2III(m 7T ,m 7T ,m K ) 
+36\/2III(m 7T , itlk, nin) — 36v / 2J//(m vr , m v , itlk) 
+72V2III(niK, rn v , m n ) + 48V2III(m v , itlk, m w ) 
+48v / 2///(m 7r , rriK, rn v ) + 24V2III(niK, fn n , m^) 

-12v / 2777(m r ,,m 7r ,mx)]. (59) 
Form Figure 4, we obtain 

M A ,K,3- P oint = j^eG F Fe v sm9{6\^m K A(rn K ,rn K ,m K ) 

+V2(m 2 K + ml)A(m n , m n , m K ) + ^-(m\ + m\)A(m K , m n , m n ) 
^~-(3m\ - ml)A(m K , m„ m v ) + ^{ m \ - ml)A(m n , m v , m K ) 
~-^( m x - ml)A(m K , m v , m n ) 

-V2(m 2 K + m^K, m n , m K ) - 6V2m K B(m K , m K , m K ) 
-^Y^k + m 2 n )B(m n , m K , m n ) - ^-(3m\ - ml)B(m v , m K , m v ) 

+ ^r( m K ~ ml)B(m n , m K , m v ) - ^^(m 2 K - m\)B(m K , m n , m v ) 
\J~2 \/2 

+ -^-( m K - m l) B ( m V' m K, ™7r) + ^~( m X ~ m l) B ( m V m K) 
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+V2C(m K , m K , m K ) —C(m w , m,, m K ) —C(m K , m T , m^) 

h\J~2 7\[2 
-7V2D(m K ,m K ,m K ) —D(m«,m«,m K ) —D(m w ,m K ,m n ) 

V2 V2 r 
— —D(m v ,m n ,m K ) —D(m n ,m K ,m v ) - 2\ / 2D(m K ,m n ,m^) 

— —D(m v , m K , m n ) —D(m v , m K , m^) 

-7V2E(m K ,m K ,m K ) —E(m 1T) m K ,m 1T ) —E(m K , m n , m n ) 

— — ^(m^m^mj - —E(m K ,m v ,m 7T ) —E(m K , m v , m v ) 

— —E(m K , m w , m v ) - 2V2E(m n , m v , m K ) 

\[2 \[2 
+V2F(m K ,m K ,m K ) - —F(m 7 , ) m K) m 1T ) - —F(m 7 ,,m 1T ,m K ) 

3^/2 \[2 
—F(m v , m v , m K ) - —F(m v , m K , m n ) - V2F(m v , m n , m K ) 



\[2 

-V2G(m K , m K , m K ) - V2G(m n , m K , m n ) + —G(m K , m n , m n ) 

3v / 2„ / x V2„, s V2„, x 

H — —G{m K , m v , m v ) + —G(m K , m v , m n ) + — G(m K , m n , m v ) 

H — — G(m n , m n , m K ) + V2G(m n: m K , m v ) - —G(m^ m v , m K ) 

-V2H(m K ,mK,'m K ) H — H(m n , m n , m K ) - V2H(m K ,m 7T ,m n ) 

V2 r V2 

+—H(m n , m K , m n ) + V2H(m K , m n , m v ) + —H^, m K , m v ) 

3^/2 V2 V2 
H — —H(m v , m K , m v ) + —H(m v , m K , m n ) —H(m v , m n , m K ) 



+8V2I(m K , m K , m K ) + \/2I(m K , m n , m,) + \p2l[m 1 ,, m K , m n ) 
H — —I{m n , m n , m K ) —I(m v , m,, m K ) + V2I(m K , m v , m n ) 

r V2 V2 
+V2I(m T) , m K , m n ) —I(m n , m v , m K ) -^I{ m m rn n , m K ) 
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-V2J(m K , m K , m K ) - V2J(m n , m n , m K ) + —J(m K , m„ m w ) 
— 2~ Jym^.m K ,m v ) + —^-J(m n ,m K ,m n ) + — J(m K , m^m^) 
+—J(m K ,m v ,m n ) + -^-J{m K ,m m m^) + \/2J(m^,m v ,m K ) 

+8\/2K(m K , ttik, m K ) + V2K(m K , m n , m n ) + V2K(m n , m n , m K ) 
H — ^ Kym-x, m K , —K(m v , m K , m v ) - —K(m w , m K , m v ) 

— —K(m v , m K , m n ) + V2K(m v , m w , m K ) + V2K(m K , m w , m v ) 

-\/2L(m K , m K , m K ) - \/2L(m K , m n , m,) + —L(m v , m w , m K ) 

H — — -L(m^, m,,, m^) + —L(m v , m n , m K ) + V2L(m K , m v , m^) 

3 \/2 \/2 
H — —L{m^ m K , m n ) + —L(m n , m v , m K ) —L(m v , m K , m^)} (60) 



is/2 

-eG F Fe u cos 9{V2(m 2 K + ml)A(m K , m n , m K ) 



12F 3 

20\/2 2 2v^2, 2 2\ * / \ 

H — m^fm^m^mj + —^(m K + mjA(m n ,m K ,m K ) 

%/2 

— - m 2 K )A(m K , m K , m v ) 

-\p2{m 2 K + m 2 n )B(m K , m K , m n ) — (m 2 K + m\)B[m K , m n , m K ) 

20y / 2 2 2y / 2 2 2 



m^B{m^, m n , m n ) + —-(m^ - m K )B(m v , m K , m K ) 



^ TV' V " ' " ' 71 / 1 ^ 

-^(ml - m 2 K )B(m K , m K , m v ) 

+V2C(m K ,m 7T ,m K ) - ^^C(m n ,m 7T ,m 7T ) - < ^^C(m n ,m K ,m K ) 
' '^C(m K , mA-, m*) - ^C(m K , m K , m v ) 



2 v ' ' ' 2 



-D(m K , m K , mj —D^m^, m n , 



2 



1 

-D(m K , rn n , m K ) —D(m K , %, m^) - 2V2D(m v , m K , m K ) 
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■y%if, rn n , m K ) — E(m n , m n , 

10V2 



r y/2 
—E(m n , m K , m K ) - 2V2E(m K , m K , m v ) —E(m K , m v , m K ) 



-—F(m n , m K , m K ) — F(m 7r , m w , — F(m K , m K , m w ) 



-V2G(m K , m n , m K ) + — ^— G(m^, m^, m^) + — — Gfm^, m x , m x ) 

H — —G(m K , m K , m,) —G(m K , m K , m v ) + V2G(m K , m v , m K ) 

2\/2 4\/2 
-V2H(m n , m K , m K ) + -^-H(m K , m n , m K ) + — fffm,, m„ mj 

H — —H(m K , m K , m n ) + V2H(m v , m K , m K ) - —H(m K , m K , m v ) 

+v / 2/(m 7r , m K , m K ) + V2I(m K , m„ m# ) + : ^^/(m 7r , m^, m^) 
H — —I(m K , m K , m n ) + V2I(m K , m v , m K ) + V2I(m v , m K , m K ) 

4:\/Qi 2\/2 

-V2J (m K , m K , m n ) + -^-J(m n , m n , m n ) + -^-J(m n , m K , m K ) 

h — y~ J ( mx ' mn > mK ^> + ^ J ( mK > mK > m »?) — Y J ( mK ' mr " mK ^ 

+V2K(m K , m K , m n ) + V2K(m n , m K , m K ) + K(m n , m n , m n ) 

H — —K(m K , m n , m K ) + V2K(m v , m K , m K ) + V2K(m K , m K , m v ) 
o 

4"\/2 2\/2 
-v / 2L(m 7r , m^, m^) + — ^— L(m 7r , m^, mj + —^-L(m K , m K , mj 

3v^2 a/2 
H — —L(m K , m n , m K ) —L(m K , m v , m K ) + V2L(m v , m K , m K )\Ql) 
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Appendix B 

We list the functions {II, III, • • • , A, B, • ■ • , L} in Figures 3 and 4 by scalar integrals of 
Paia2a3 anci one-loop tadpole integrals T\, T 2 in the Euclidian space. We note that the 
function / does not contain g^ v and thus it has no contribution to Fa- For simplicity, we 
only give the formulas related to terms with due to the definition of Fa- We have 



II(mi,m 2 ,m 3 ) 



(2n) 2n (n-l) 
-p p iii( m i, m 2,rn 3 ) 



Ti(m 2 )Ti(m 3 ) - m 1 P ul (m 1 ,m 2 ,m 3 ) 



+ 



(62) 



III(m 1 ,m 2 ,m 3 ) 



-ig 



(2n) 2n (n- 1) 



1 



1. 



T^m^TxK) - -T^m^mi) - -Ti(m 2 )Ti(m 2 ) 



+ ( m i + m 2 -m 3 - £ )P lll (m 1 , m 2 , m 3 ) - P m (m,i, m 2 , m 3 ) 



A(mi,m 2 ,m 3 ) = 



-PJW (mi , m 2 , m 3 ) - (mi , m 2 , m 3 ) 



£ 2 " 



+ 



(63) 



-2<y 



/ii/ 



(27r) 2 «(n - 1) 

-^2- P 2 2 i°i(^l>^2,m 3 ) 



| [Pi°i°i(mi , m 2 , m 3 ) - miP 2 i°i(mi, m 2 , m 3 ) 



+ 



(64) 



B(m 1 ,m 2 ,m 3 ) 



(2*) 



2((/ 



——^ J dx(-l) [T 2 (m2)Ti(m 2 ) -T 2 (m2)Ti(m2) 



-P m (m 2 ,mi,m 3 ) + (m x + m 2 - m 3 - £' )P 211 (m 2 ,mi,m 3 ) 
-2P 2 1 i° L (m 2 , mi, m 3 ) - 2P 2 ° 1 1 1 (m 2 , mi, m 3 ) 

r P 2 ii(m 2 ,mi,m 3 ) 



f 2 ' 



+ 



(65) 



C(m 1 ,m 2 , m 3 ) + P(m 3 , mi, m 2 ) 



(2tt; 



2ri , 



n 



1) 



2(1 - x)(p M ■ k)P™ 1 (m 1 ,m 2 ,m 3 ) 



Pi\°iK,m 2 ,m 3 ) - Ti(m 2 )Ti(m 2 ) + Ti(m 2 )Ti(m 2 ) - T^mfiT^mg) 



£' 2 



+2P 1 1 1 ° 1 (m 1 , m 2 , m 3 ) + (-m 2 + m 2 + m 2 + £ ,2 )P 1 ° 1 ° 1 (mi, m 2 , m 3 ) 

(v ■ £') 

-2m 2 (l - x)(p M ■ k)P™ 1 (m 1 ,m 2 ,m 3 ) + 2m 2 ^|^ — P 2 l ° 1 (m 1 ,m 2 ,m 3 ) 

+m 2 T 2 (m 2 )Ti(m 2 ) - m 2 Pi°i°i(mi, m 2 , m 3 ) + m 2 T 2 (m 2 )Ti(m 3 ) 
-2m 2 P 2 1 1 ° 1 (mi, m 2 , m 3 ) - m 2 (-m 2 + m 2 + m 2 + £ /2 )P 2 ° 1 ° 1 (mi, m 2 , m 3 ) 
2(1 - x){pm • k) 20 2(p M • f ) 30 



£/2 



-P 211 (mi,m 2 ,m 3 ) + 



f 4 



P 211 (mi,m 2 ,m 3 ) 



2 1 

-— P 2 3 n( m i 5 m 2, - ^( _m i + m 2 + m 3 + ^' 2 )P 2 ii(m 1 ,m 2 , m 3 ) 



1 

f 2 ' 



+ ^1(^2) 



d n p(p 



ff'\2 



1 



(p 2 + m 2 ) 2 f 2 



— P m (mi,m 2 ,m 3 ) 
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'A2 



+ 



(66) 



D(m 2 , mi,m 3 ) + E(m 1 ,m 3 , m 2 ) 



(2vr) 2 "(n - 1) 



dxlT^mj) / d n p 



■pr- 



ill 



(p 2 + m 2 ) 2 

Pi i\(m 2 , mi, m 3 ) + (1 - • fc)*m(™2, mi, m 3 ) 



-(l-x)(f • A;)T 2 (m 2 )Ti(m 2 ) 



(1 - x)(f • k) 01 . , (l-a:)(^-fc)„io/ 

-^Oms, mi, m 3 ) + ^ -^mK, mi, m 3 ) 



f 2 



^(m 2 )^™ 2 ) + ^mlP^m^m^mz) - ^(m^T^m 2 ) + ^Ti(m 2 )Ti(m 2 ) 



--(m 2 - £' 2 - m 2 3 )PZ(m 2 , mi, m 3 ) + P^m,, mi, m 3 ) + (p M • 0^2 K^m 2 ) 



+ (m 1 +m 2 — — m 

(1 - x)(£' ■ k) 
F 2 



(pm ■ n p o li(m2? mij mg) _ (1 _ . fc)P oo i(m2) TOi> mg) 

P 2 ii(m 2 ,mi,m 3 ) - — — — -P 2 1 1 ° 1 (m2,mi,m 3 ) 



f 2 



I' 2 



+ ip i ° 1 1 (m 2 ,mi,m 3 ) - ^^"(mj, mi, m 3 ) + ^(m^T^m 2 ) - iT 2 (m 2 )Ti(m 2 ) 
+i(m 2 - £' 2 - m^P^Vms, m 1 ,m 3 ) - P^^, m 1 , m 3 ) 



(Pm • p o2 i(m2; ms) + 2(1 _ ^ . fc ) P oi i ( m2> mi5 ms ) 

P 2 02 i(m 2 , mi, m 3 ) + — — ^— - ^ P 2 n(m 2 , mi, m 3 ) 



f 2 

2(1 - • A;) 



P^K, mi, m 3 ) + m\P*\ x {m 2 , mi, m 3 ) - £ ,2 T 2 (m 2 )Ti(m 2 ) 
(m 2 - f 2 - m^P^ms, mi, m 3 ) + 2P 2 02 1 (m 2 , mi, m 3 ) 
o(PM-f) 



f 2 



-P 2 ii(m 2 , mi, m 3 ) + 2(1 - x)(f • A;)P 2 1 1 ° 1 (m 2 , m 1; m 3 ) 



2(1 - x)(f • A;) n , 2(1 - x)(£' ■ k) 20 
+- Lp 2n(™ 2 , mi, m 3 ) + -5= -P 2 Ti(™2, mi, m 3 ) 

-Pii°i( m 2, mi, m 3 ) + m\P\l x (m 2 , m 1 ,m 3 ) - (ml - £' 2 - m^P^m^ mi, m 3 ) 



+2P 21 1 1 (m 2 ,mi,m 3 ) - 
, 2(1 - x)(f • A;) Dll 



P 211 (m 2 ,mi,m 3 ) 



f 2 

2(1 - • &;) 

f 4 



f 4 
•?ii 

P 2 2 ii("72,mi,m 3 ) - ■^P 1 1 1 1 1 (m 2 ,mi,m 3 ) 



P 2 n(m 2 ,mi,m 3 ) + 2 



(l-x)(f .fc) Dl2 , . 

^P 2 Yi(m 2 ,mi,m3) 



f 4 



2 7~»1 1 / \ ( 1 ^ /?? o ) ii , > 

+"^2 m 2 p 2ii( m 2, mi, m 3 ) ^ P 2 n( m 2, mi, m 3 ) 

2 1 

+^A 1 i 2 iK,mi,m3) - ^Ti(m^) 



/\2 



• d"p(p • £') 
(p 2 + m 2 ) 2 



+ 



(67) 



H(m 2 , mi, m 3 ) + J(m 1 ,m 3 , m 2 ) 
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-9T fj [2(p M .t') 



(27r) 2 ™(n- 1) 



/ dx | VM (ii2 P m( m 2- m u m s) + 2^11! (?ri2, mi, m 3 ) 

-2(1 - x)(t ■ fc)Pg , 1 (m 2 ,mi,m 3 ) - ^"^'^ ^iKmi.m,) 
2(1 - a;)(f • fc) 01 _ _ . , m( ,. rr( ,. 



■P m (m2,mi,m 3 ) + 7i (m^T^m^ 



£/2 

+(-m 2 + m 2 - £' 2 - m^P^m^ m u m 3 ) - 2P 1 1 1 ° 1 (m 2 , mi,m 3 ) 
-2<(m2,mi,m 3 ) + 7\(m 2 )7\(m 2 ) - Ti(m 2 )7\(m 2 ) 
2m2(p M .f) D0 o_2 D io 



f 2 



-P 211 (m 2 ,m 1 ,m 3 ) - 2m 2 P 211 (m 2 ,m 1 ,m 3 ) 



+2ml(l - x)(£' ■ k)PZ(m 2 ,m u m 3 ) + ^K 1 W ' k ) pio^ ^ 
+ 2m 2 (l -x)(£> ■ k) p o li{m2 mi mz) _ m 2 Ti{m 2 )Ti{m 2 ) 

-m\{-m\ + m\- £' 2 - ml)P^ 1 (m 2 , m 1 ,m 3 ) + 2m 2 2 P 2 1 ° 1 (m 2) m 1 , m 3 ) 
+2m\P 2 l 1 (m 2 , m 1 ,m 3 ) + m 2 2 T 2 {rr^ 2 )Tx{m 2 3 ) - mlP^m^ mi,m 3 ) 
2(Pm • f ) d21 , , 2 d30 



P 2U (m 2 ,m 1 ,m 3 ) - —P 2ll {m 2 ,mi,m 3 ) 



£,4 ~ ZLL\""*1 ""D ■■"*) £i 2 

+ 2(1 - xW ■ k) p| o (w2> TOi> + 2 (1 - x)(f • fc) p 3o (m2; Wi> ma) 



£/2 



(p 2 + m 2 ) 2 



;-m| + m\ - £ 12 - m 3j 20 ... , 



P 211 (m 2 ,mi,m 3 ) + -^P 211 (m 2 ,m u m 3 ) 

« P ^75 



£/2 -^llv z ' J-' °' 1 £/2" 



21 r 

-— P 2 2 1 1 1 (m 2 ,mi,m 3 ) - j^Piu(m 2 ,m 1 ,m 3 ) + 7\(m 2 ) y ■ 



(p 2 + m 2 ) 2 
+ ••• , 



J(m 3 , m 2 , mi) + K(m 2 , m u m 3 ) 

(2^(1-1) / ^ • ^(mlWm!) 

-(pm • l')Pii 1 {m 3n m 2 ,mi) - VM ^ 2 Pi\°i(m 3 , m 2 , mi) 



f2 P m (m 3 ,m 2 ,mi) + — P^m^m^m^ 

-ip^m^TxK) + ^(m^TiK) + ^(m^T^m 2 ) 
~^( m 2 - m i - + m 3)^i°Li( m 3, m 2 , mi) + Pm^, m 2 , mi) 



+P 1 1 1 1 (m 3 ,m 2 ,mi)-(l-a;)(£'.A;)T 2 (m 2 )Ti(m 2 )+Ti(m 2 ) / 
+ (m 2 + m 2 - £' 2 - m 2 ) x [(p M • 0*m(™3, ^2, 

^ -4°i( m 3, m 2, »i) + ^ M m o^ P 2 li{m^ m 2 , mi) 



£/ 2 211V "3> "1/ 1 £/ 2 
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(I — x )(£' . fa) | 

^P^ 1 (m 3 ,m 2 ,m 1 ) + -T 2 (m 2 )Ti(m 2 ) - -P^mg, m 2 , mi) 



f 2 



1 

2" 



|r 2 (m 2 )Ti(m 2 ) + 2 : (m 2 + m 2 - £' 2 - m^P™^, m 2 , rrn) 



-4°i( m 3. m 2. m l) 



P 2 1 1 1 (m 3 ,m 2 ,m 1 )] 



2(p M ■ £') 

-2(p M ■ £')P 2 u(m 3 , m 2 , mi) — P^m^ ™2, mi) 

2(pM-f) 



£/2 



P 2 ii(m 3 ,m 2 ,mi) + 2 



1 



-P 2 1 1 1 1 (m 3 ,m 2 ,mi) 



£/2 -"JA °' z ' ^ 1 £/2 

+Pii°i(m 3 , m 2 , mi) - (m 3 + m 2 - £' 2 - mi)P 2 i°i(m 3 , m 2 , mi) + 2P 2 2 1 ° 1 (m 3 , m 2 , mi) 

2fo • £0 

+2P 2 1 1 1 1 (m 3 , m 2 , mi) - 2(p M • O-^ii^s, ™2, mi) ^ ^nK, m 2 , mi) 



2(pm 



P 2 n(m 3 ,m 2 ,mi) + 2 



£/ 2 211V "Oj »z, ly 1 ^ /2 

-f 2 T 2 K)Ti(m 2 ) - (m 2 3 + m 2 - f 2 - m^P^m,, m 2 , mi) 



P 2n (m 3 , m 2 , mi) + P in (m 3 , m 2 , mi) 



-2 
+2 



ui , o 002 1 . 2(p M • (!) nl , 

211 

-P 2 2 ii( m 3, Ws, ™i) - ^ ^./^ anK, ™2, mi) 



+2P 21 1 1 (m 3 , m 2 , mi) + 2P 2 02 1 (m 3 , m 2 ,m 1 ) 



f 2 



P 2 Yi(m 3 ,m 2 ,mi) 



f 4 



— P 21 2 i(m 3 , m 2 , mi) + ^Pii\(m 3 , m 2 , mi) 



f 4 

(m 2 + m 2 



f 4 
1 

^2"- 



f 2 



m 2 ) 2 
— L ^2iiKi m 2 , mi) + —Pll^m^ m 2 ,m x ) 



+ 



(" 2 



P 2 1 ii(m 3 ,m 2 ,mi) - ^Ti(m 2 ) ^ 



/\2 ' 



(p 2 + m 2 ) 2 



+ 



(69) 



G(mi, m 2 , m 3 ) + L(m 3 , m b m 2 ) 



I dx |(pM • A;)(l - x)T 2 (m 2 )Ti(m 2 ) - T 2 (m 2 ) | 
-Pii°i(mi, m 2, m 3 ) - (p M -k)(l- x)P^ 1 (mi,m 2 , m 3 ) 



d n pp 2 



(2n) 2n {n- 1) 



(p 2 + m 2 ) 



f 2 



+P 1 ° 1 1 1 (m 1 ,m 2 ,m 3 ) + iTi(m 2 )Ti(m 2 ) - ^(m^T^m 2 ) + ^(m^T^m 2 ) 



+ 2 ( m i ~~ m l ~ ^ ~ m l)Pni( m i, m 2, m 3 ) - PiW^i, ™2, m 3 ) 
-Pi i\(mi, m 2 , m 3 ) - (1 - x)(p M ■ £;)T 2 (m 2 )Ti(m 2 ) - (m 2 + m 2 - £' 2 - m 2 ) x 
(Pm ' 

-4°i( m i) m 2, m 3 ) - (p M • - ^)-P 2 i°i( m i> m 2, m 3 ) 



f 2 



+P 2 n(mi,m 2 ,m 3 ) + -T 2 (mi)Ti(m 2 ) - -P m (mi, m 2 , m 3 ) + -T 2 (mi)Ti(m 3 ) 
+ 2~( m i ~ m 2 - - m t) P 2ii( m i, m 2, m 3 ) - Pan (mi, m 2 , m 3 ) 



-P 2 u 1 1 1 (mi,m 2 ,m 3 )]+2 



f 2 



P 211 (mi,m 2 ,m 3 ) 



-2(p M • fc)(l - x)P 2 1 1 ° 1 (mi, m 2 , m 3 ) + 2P 2 1 1 \(m 1 , m 2 , m 3 ) 
-Pii°i(mi, m 2 , m 3 ) + (m 2 - m 2 - f 2 - m^P^mi, m 2 , m 3 ) 
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-2P 2 2 1 ° 1 (mi, m 2 , m 3 ) - 2P 2 \ 1 1 (mi, m 2 , m 3 ) 

+2 PM ^ Pmipii, m 2 , m 3 ) - 2(p M ■ k)(l - x)P^ 1 (m l , m 2 , m 3 ) 
+2P 2 02 1 (m 1 , m 2 , m 3 ) - P^K, m 2 , m 3 ) - ^(m^T^m 2 ) 

+ (m\ - m 2 - £' 2 - m^P^TOi, m 2 , m 3 ) - 2P 2 1 1 \(m 1 , m 2 , m 3 ) - 2P^ 1 (m 1 , m 2 , m 3 ) 

+ ^jPM_n P 2i i(mi) m2; ms) _ 2 ^'^-^ i (mii TO2> ros) 

2 1 

+ -^P 2 n( m i, m 2,m 3 ) - —P^ 1 (m 1 ,m 2 ,m 3 ) 

(m? — m 2 — £' 2 — mo) n 2 21 

+ ^2 P 2 n{mi,m 2 ,m 3 ) - —P 2n (m u m 2 ,m 3 ) 



2 „12 / x 1 ^ / 2x f d n p(p-£') 2 I 

■^P 211 (m 1 ,m 2 ,m 3 )- ; ^T 1 (m3) j ^ - ^ 2)2 j> + ■ ■ ■ , (70) 



where = £ + kx, n = 4 — 2e, Pm represents the meson momentum, and {• ■ •} correspond 
to the terms without g^ v . 
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Appendix C 



In this Appendix, we will replace the set of ten functions P 2 < f 1 (m 1 ,m2,m3;! 2 ) by the 
following equivalent one of Hi{mi, m 2 , m 3 ; £ 2 ) [17]. The functions Hi are free of quadratic 
divergencies and for this reason they have simpler integral representations. 

For the well known one-loop integrals with the definition 7 S = 7 — 1 — ln(47r), we have 



ji 2\ _ 4-D f d D q 



m 



m 



IQlT 2 \4lTfJ, 2 



(2n) D q 2 - m 2 

r(-i + e) 



-m 2 II / \ 

i6^U- 7s - ln \T 2 j +e 



7T 

12 ~ 7s 



, / m 

1+ln U 



,(71) 



(27r) 4 -"T 1 (m 2 ) = (27r/i) 4 -™ / d n V -^- 

J p A + 



2 2 , 
-m 7r < - 



7s 



.2 
12 



7T 

T7T -7s 



ln bJ + 2( 7s+1+ln b 



,(72) 



(27r) 4 -"T 2 (m 2 ) = (2^)^™ / tTp— ^ 



(p 2 + m 2 )' 



m 



1 /m 2 \ , /m 2 * 

+ 5 ln b xln b. 



7T 1 it , I m \ , ( m z 

_ + _ +7 , + | + ln /_j +7sln f_ 



(73) 



We note that the expressions of Eq. (71) and Eqs. (72) and (73) are defined in Minkowski 
and Euclidian spaces, respectively. 

For the two-point functions, the relations between Pffi and Hi are given by 



P^ 1 (m 1 ,m 2 ,m 3 ; 



l —{{m\ + f)H 1 (m 1 ,m 2 ,m 3 ) + H 2 (m 1 ,m 2 ,m 3 ) 

n — 6 

+m 2 2 H 1 (m 2 , mi ,m 3 ) + m^K, m u m 2 )} , 



(74) 



P^ 1 (m 1 ,m 2 ,m 3 ;£' 2 



lll (m 2 ,m 1 ,m 3 ;£ 



pOl 

-1 
n — 



-77^111(^1,^2,^3) - P 2 ii(m u m 2 ,m 3 ) 



-mlH 2 {m u m 2 , m 3 ) - ^(m, m 2 , m 3 ) 
-m 2 3 H 3 (m 3) m 2 ,m l ) , 



(75) 
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P 111 (m 1 ,m 2 ,m 3 ;£ ) 



-7T [^2^211(^2, m 1 ,m 3 ) - m 2 3 P 2 ^(m 3 , m 2 , mi) 

72 — 2 L 

- m 3A n( m 3, m 2 , mi) - £ 2 777 2 P 2 ii(m 3 , m 2 , mi) 
+m^P 2 1 1 1 1 (m 1 ,m 2 , m 3 ) - P£ 1 (m 1 ,m 2 , m 3 ) 



P^^m^ma^g) 



(76) 



P lu (m 1 ,m 2 ,m 3 ;£ ) = P m (m 2 , mi, m 3 ; £ ) 



m 2 P 2 2 i°i( m 2, mi, m 3 ) + 772^2^(7723, m 2 , mi) 



02 



+m 1 P 211 (mi,m 2 ,m 3 ) - P 2U (m l ,m 2 ,m 3 ) 



(77) 



For the three-point functions, we use [17] 



P 2 ii(mi,m 2 ,m 3 ; £ ) = H 1 (m 1 ,m 2 ,m 3 ) , 



(78) 



P 211 (m 1 ,m 2 ,m 3 ; 



. 0^ 



-H 2 (m 1 ,m 2 ,m 3 ) - £ 2 H l (m l ,m 2 ,m 3 ) , 



(79) 



P 2 ° 1 1 1 (mi,m 2 ,m 3 ;. 



-H 3 (m 1 ,m 2 ,m 3 ) , 



(80) 



* r r 1 

P 2 ° 1 (m 1 ,m 2 ,m 3 ;£ 2 ) = H 4 (m 1 ,m 2 ,m 3 ) H { (n - l)£ 2 - m\ H 1 (m 1 ,m 2 ,m 3 ) 

n L L J 

+2(72- l)if 2 (mi,m 2 ,m 3 ) +P 1 1 ° 1 (mi,m 2 ,m 3 )} , 



(81) 



P^ 1 (m 1 ,m 2 ,m 3 ;£ 2 ) 



H 5 (m 1 , m 2 , m 3 ) + £ 2 H 3 {m x , m 2 , m 3 ) 
I 2 

+ — (m 2 + m 2 - m 2 + £ 2 )H 1 (m 1 ,m 2 ,m 3 ) 

+2H 2 (m 1 ,m 2 , m 3 ) - P™ 1 (m 1 ,m 2 , m 3 ) 
+T 2 (m 2 )Ti(m 2 ) -T 2 (m 2 )Ti(m 2 )} , 



(82) 



P 211 (mi,m 2 ,m 3 ;£ ) = H 6 (rrii, m 2 , m 3 ) H -m 2 H 1 (m 1 ,m 2 ,m 3 ) 



n 



+T 2 {m 2 l )T l {ml) 



(83) 



P|i°i(mi,m 2 ,m 3 ;r) = 



-H 7 {m 1) m 2) m 3 ) - P^ 1 (m 1 ,m 2 ,m 3 ) 
3£ 2 (£n-l 



77 + 2 



£ 2 -mn^(mi,m 2 ,7n 3 ) 



+ 



n - 1)£ 2 - m 2 i7 2 (mi,m 2 ,m 3 ) + nP 4 (mi,777 2 ,m 3 )| (84) 



P 2U (m 1 ,m 2 ,m 3 ;£ ) = -H s (m u m 2 ,m 3 ) - 



3£ 2 
n + 2 



(n - l)P 5 (777i,7?7 2 ,777 3 ) 
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+ 



n — 1 
3 

n — 1 



£ 2 -m{) H 3 (m 1 , m 2 , m 3 ) - ^(mi, m 2 , m 3 ) 



; m 2 + m 2 _ m 2 + p) Hl ( mu m2? m3 ) 



n(n + 2) 

+2H 2 (m 1 ,m 2 , m 3 ) - P 1 c 1 1 (mi,m 2 , m 3 ) 
+T 2 (m?)Ti(7r^) -T 2 (m?)Ti(ml " 



(85) 



P 2U (m 1 ,m 2 ,m 3 ;£ ) = -Hg(m 1 ,m 2 ,m 3 ) - £ H 6 (m 1 ,m 2 ,m 3 ) 

'2£ 2 



2H 5 (m 1 ,m 2 , m 3 ) + I — - m 2 j H 2 (m 1 , m 2 , m 3 ) 



n + 2 

+ (m? + r^-ml + ^ 3 (m 1 ,m 2 ,ms) 
+P 1 ° 1 1 1 (mi,m 2 ,m 3 ) 



{-Pi i°i(^i,^2,m 3 ) 



n(n + 

m 2 _ (n + 1)m 2 _ m 2 + £ 2] 

+T 2 K)T 1 (m2) - (n - l^m 2 )!^ 2 )} , 



(86) 



P 2n (mi,m 2 ,m 3 ;£ ) = -Hi (mi,m 2 ,m 3 ) 

+e 2 T 2 (mj)T 1 (m 2 3 ) 
The functions of i/j are expressed as follows 



3£ 2 



n + 2 



-m\H 3 {m x ,m 2 ,m 3 ) 



(87) 



H x {m x ,m 2 ,m 3 ;£ 2 ) = rr 4 



2 1 1 7T 2 

^2 ~~ A 7mi / ~~ 2 ^ ~\2 ~ 7 
+/ii(mi,m 2 ,m 3 )] , 



mi 



H 2 (mi,m 2 ,m 3 ;£ 2 ) = n 



4«2 



2 1 J _ , 13 7T 2 7mi 2 

-77" -7, 



— + -(- - 2 7mi ) + — - — + 



A 2 A v 2 
-/i 2 (mi,m 2 ,m 3 )] , 



8 12 



mi 



(89) 



if 3 (mi,m 2 ,m 3 ;£ 2 ) = 7r 



4/)2 



1 11 13 f 2 

A 2 ~ AM 7mi) ~ i~6 + 24 
+ /i 3 (mi,m 2 ,m 3 )] , 



7m 1 . 7 



mi 



(90) 



H 4 (m 1 ,m 2 ,m 3 ;- 



4/)4 



7T t 



3 1 3 7mi 



2A 2 A 2 



175 vr 2 3 7 , 

1 h 

96 16 



2 

mi 



+ -^h 4 (mi,m 2 ,m 3 ) 



H 5 (mi,m 2 ,m 3 ;£ 2 ) = 7r 4 £ 4 



13 7mi , 175 vr 2 3 7 ^ 



+ 



4A 2 A 4 



192 32 



(91) 



(92) 
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H 6 (m 1 ,m 2 ,m 3 ;£ 2 ) = ir 4 £ 4 



H 7 (m 1 ,m 2 ,m 3 ;f) = n 4 £ 6 



2A 2 A l 24 2 ' 32 48 24 



7 3 
+ + ^e( m i, m 2 , m 3 ) 



J_ x 287 7T 2 5 7mi 

'A 2 A l 24 7mJ 192 24 24 

7 2 1 

-y~ - -h 7 {m 1 ,m 2 ,rn z ) 



H 8 (mi,m 2 ,m 3 ;£ 2 ) = ir A f 



1,5 7 m 287 7T 2 57 



H — 7~r + 



mi 



H 9 (mi,m 2 ,m 3 ;£ 2 ) = n 



4/)6 



— + — (— + — ) 

2A 2 A v 48 2 ; 384 ' 48 ' 48 



7 1 
+ ~f- + 2^8(^1,^2,^3) 



— ( — + l^L) 95 7 mi 

v 1 " 3 j 192 72 24 



Hi (mi,m 2 ,m 3 ;£' 2 ) = 7r 4 £ 6 



3 A 2 A v 24 
7 2 1 

-y- - -/i9(mi,m 2 ,m 3 ) 



1 , 1 , 1 + Tmi, _ 283 + 7T 2 + 7mi 



4A 2 A v 96 4 7 768 ' 96 ' 96 

7 2 1 
+ "f 1 + 2^10(^1,^2,^3) 



where A = — 2e and 7 m = 7 + ln(7rm 2 //x 2 ). 

The ultraviolet finite parts hi(mi, m 2 , m 3 ) of the function H i (mi,m 2 , m 3 ; £ 2 ) 
following one-dimensional integral representations: 



h 1 (m 1 ,m 2 ,m 3 ) 


= / dy[g(y)], 
Jo 




h 2 (m 1 ,m 2 ,m 3 ) 


= I' dy[g(y) + A(y)] , 
Jo 




h 3 (m 1 ,m 2 ,m 3 ) 


= f dy[g{y) + A(y)](l - y) , 
Jo 




h 4 (m 1 ,m 2 ,m 3 ) 


= I' dy[g(y) + A(y) + f 2 (y)] , 
Jo 




h 5 (m 1 ,m 2 ,m 3 ) 


= C dy[g(y) + h{y) + / 2 (y)](l 
Jo 


-y), 


ha(m 1 ,m 2 ,m 3 ) 


= C dy[g(y) + h{y) + / 2 (y)](l 

JO 




/t 7 (mi,m2,m 3 ) 


= ^ dy[g(y) + A(y) + f 2 (y) + 
Jo 


km , 
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h 8 (m 1 ,m 2 ,m 3 ) = f dy[g(y) + + f 2 (y) + f 3 (y)](l - y) , 
Jo 

h 9 (m 1 ,m 2 ,m 3 ) = f dy[g(y) + f^y) + f 2 (y) + / 3 (y)](l - yf , 
Jo 

/i 1 o(m 1 ,m 2 ,m 3 ) = f dy[g(y) + /i(y) + f 2 (y) + f 3 (y)](l - yf . 
All ten integral representations are built up by the following four basic functions: 



(98) 



g{y) 
h(y) 
/ 2 (y) 
/s(y) 



^ 1^—1 + ^ f — I + 2/1 ln I -^-r I + 



1 
2 

1 

3 

1 

4 



.1-3/1, 

l-^ 2 

/t 2 



+ y 2 ln 



. 1 - 2/2 y 
f 2/i 



1/2 ln 



/ 1/2 



1/2 In 



1/2 



V 1/2 - 1/ ' 



2 1-z/ 2 /l-z/ 2 \ 2 



/t 2 2/t 2 



1/1-1/ Vl/2- 1 

+ l/i 3 In 



3i„ / 2/1 \ , „. 3 ,_ / 2/2 



2/1-1 



+ 1/2 ln 



,2/2-1 



4 /l 3 \ 1 - v 2 1 /l - z/ 2X 2 
n - - + 



k 2 V3 k 2 / 2/t 2 2 V k 2 



K 



2\ 3 



where 



+ l/i In 



2/1,2 = 



z/ = 



4i / 2/1 \ , 41../ 2/2 



,2/i-l 



+ l/ 2 4 ln 



,2/2 - 1, 



(99) 



ln(l - i) 



dt, 



1 + /c 2 - v 2 ± ^(1 + k 2 - z/ 2 ) 2 + 4z/ 2 k 2 - 4m 2 ry 



2/t 2 



ay + 6(1 - y) 



m 2 t m 3 2 
' a = ~7 2 ' = 2 ' K = 2 ' 



(100) 



y(l - y) m\ m\ ml 

Finally, we must transform the parameters back into the Minkowski space and change 
the inward directions of the momenta /, k for the final particles by outward ones: 

f - -(m* -2p-k), 



e -> -[m 2 -2p-k(l-x)}, 



p ■ k — > p ■ k , 



p-H -> m 2 K ^-p-k(l-x) , 



k ■ £' -> -p-fc, 



(101) 
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where < p ■ k < (m 2 Kn — mf)/2, and we have replaced d 2 by — d 2 as calculated in the 
Euclidian space. 
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